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Abstract 



The structures of the ground and excited states of ^^Be were studied with a 
microscopic method of antisymmetrized molecular dynamics. The theoretical 
Tlj" . results reproduce the abnormal parity of the ground state and predict vari- 

\ ous kinds of excited states. We suggest a new negative-parity band with a 

O \ well-developed clustering structure which reaches high-spin states. Focusing 

^ ' on a 2a clustering structure, we investigated structure of the ground and ex- 

cited states. We point out that molecular orbits play important roles for the 
intruder ground state and the low-lying 2%ijo states. The features of the break- 
^ , ing of a clusters were also studied with the help of data for Gamow- Teller 

^ ' transitions. 



I. INTRODUCTION 

Recently, information on the excited states of light unstable nuclei has increased rapidly 
In the excited states of light unstable nuclei, an exotic molecular structure of light 
unstable nuclei is one of the attractive subjects in experimental and theoretical research. 
For example, molecular structure has been suggested to appear in neutron-rich nuclei, such 
as ^°Be and ^^Be [|Ti^Hl^. W. Von Oertzen et al. proposed a kind of exotic clustering 
structure with a 2huj configuration in the excited states of ^^Be. However, there have been 
few microscopic studies on the excited states of ^^Be. 

Needless to say, a-cluster cores are very important in the ground and excited states of 
^^Be as well as in ^Be and ^°Be. However, we should not forget the breaking of a-clusters in 
^^Be because there are many valence neutrons around the a-clusters. The experimental data 
concerning the /3-decay strength are very useful to estimate the breaking of the 2a clustering 
structure in Be isotopes, because Gamow- Teller transitions from Li to Be are not allowed if 
the a-cluster cores in the daughter state of Be are completely ideal ones with a simple (Os)^ 
configuration. It is very interesting that the recently measured /3-decay strength from ^^Li to 
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^^Be indicates significant breaking of a-clusters in tlie excited states of ^^Be p|. Therefore, 
it is important to also study the breaking of a-clusters due to the surrounding neutrons as 
well as the development of clustering in ^^Be. 

We should point out another interesting feature in ^^Be. Abnormal parity of the ground 
state in ^^Be has been known for a long time. Namely, the spin parity of the ground state 
is 1/2+, which seems to be inconsistent with the ordinary shell-model picture in which ^^Be 
with 7 neutrons may have a l/2~ state as the ground state. It has long been a problem why 
parity inversion occurs in ^^Be. As possible reasons for this parity inversion, such effects as 
the halo structure, a core deformation, clustering and paring are suggested. For example, 
the energy gain of a positive parity state of ^^Be is discussed in Refs |[13|-[17|. However, the 



ground-state structure has not been sufficiently studied by fully microscopic calculations 
without such assumptions as the existence of core nuclei. Although the clustering struc- 
ture must be important in Be isotopes, it is difficult for mean-field approaches to describe 
developed clustering structures. With a theoretical method of antisymmetrized molecu- 
lar dynamics. Dote et al. have studied the abnormal parity ground state of ^^Be without 
assuming cores or the stability of the mean field. However, in their work [Q, the effects 
of angular-momentum projections and three-body forces were approximately estimated by 
perturbative treatments. 

Our aim is to make a systematic research on the structure of the ground and excited states 
of ^^Be based on the microscopic calculations. An important point is that the theoretical 
approach should be free from such model assumptions as stability of the mean field, and the 
existence of inert cores or clusters. First of all, traditional mean-field approaches are not 
useful to study the developed clustering structure in Be isotopes. With cluster models, the 
clustering structure of the excited states of ^Be and ^°Be has been successfully explained 
rT| , p!^JTS| -pU[ by assuming 2-a cores and surrounding neutrons. Since the assumption of 2 



a-cluster cores is not appropriate to discuss the breaking of a-cores, these cluster models are 
not sufficient for investigations of ^^Be with many valence neutrons. In fact, it is difficult 
to use them to directly calculate the experimental data of the /3-decay strength from ^^Li. 
With these models it may not be possible to describe the recently discovered excited state 
at 8.04 MeV with the strong /? transition strength from ^^Li. 

We have applied a theoretical approach of antisymmetrized molecular dynamics (AMD). 
The AMD method has already proved to be useful for studying the structures of light 



nuclei ||8|, pI)|j2T| --p4| . Within this framework, we do not need such model assumptions as 
inert cores, clusters, nor axial symmetries, because the basis wave functions of a nuclear 
system are written by Slater determinants where all centers of the Gaussian-type spatial 
part of single-particle wave functions are free parameters. In AMD studies of neutron- 



rich nuclei, we investigated the structures of Be isotopes P, |22| , |24| . In AMD calculations, 
many kinds of experimental data for nuclear structure have been reproduced. Due to the 
fiexibility of the AMD wave function, we have succeeded to describe the structure changes 
between shell-model-like states and clustering states with an increase in the neutron number. 
In previous studies P,P^ on Be isotopes, the excited states of ^^Be were not studied in 
detail because the calculations are based on the variation before a total-angular-momentum 
projection. Recently, the AMD framework has been developed to be an extended version 
based on variational calculations after a spin-parity projection(VAP), which has already 
been confirmed to be powerful for studying the excited states of light nuclei. The method 
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has been applied to a stable nucleus ^^C and to unstable nuclei: ^°Be, ^^Be and ^^Be 

iiioii- 

In the present work, the structures of the ground and excited states of ^^Be were studied 
by performing VAP calculations using the AMD method. In the next section (Sec. U), we 
explain the formulation of AMD for a nuclear structure study of the ground and excited 
states. The adopted effective interactions are briefly explained in Sec. |ITl|. In Sec.|^ we 
present the theoretical results concerning the energy levels, (3 decays and E2 transitions 
compared with the experimental data. We predict a new rotational band with a well- 
developed clustering structure which comes from 2huJ configurations. In Sec.[V|, the intrinsic 
structure and behavior of the valence neutrons are discussed. We also discuss the breaking 
of a-clusters which has an important effect on the /3-decay strength. Finally, we summarize 
our work in Sec. |V| 

II. FORMULATION 



In this section, the formulation of AMD for a nuclear structure study of the excited 
states is explained briefly. For more detailed descriptions of the AMD framework the reader 
is referred to Refs. [ITUll 



A. Wave function 

An AMD wave function of a nucleus with mass number A is a Slater determinant of 
Gaussian wave packets; 
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where the ith single-particle wave function ipi is a product of the spatial wave function (pXi , 
the intrinsic spin function x^- and the isospin function Tj. The spatial part (pXi is presented by 
a Gaussian wave packet whose center is defined by complex parameters Xu, X2i, X^i. Xii is 
the intrinsic spin function parameterized by ^j, while Tj is the isospin function which is fixed 
to be up(proton) or down(neutron) in the present calculations. Thus an AMD wave function 
is parameterized by a set of complex parameters Z = {X„j,^j} (n = 1,3 and i = 1,A). 
are the centers of Gaussians for spatial parts and the parameters ^j's determine the directions 
of intrinsic spins of the single particle wave functions. 

If we consider a parity eigenstate projected from an AMD wave function, the total wave 
function consists of two Slater determinants, 



$(Z) = (1 ±P)<I> 
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where P is a parity projection operator. In the case of a total-angular-momentum 
projection(J-projection), the wave function of a system is represented by the integral of 
the rotated states, 

$(z) = p/,^,$^AfD(z) = J dnDi;j,,{n)R{n)^AMD{z), (4) 

where the function -Dj(^^ is the well-known Wigner's D function and R{Q) stands for the 
rotation operator with Euler angle Q. 

In principal the total wave function can be a superposition of independent AMD wave 
functions. For example, a system is written by a superposition of spin-parity projected AMD 
wave functions PI^,(^amd as follows, 

$ = cP^tK'^AMDil) + dPif-^,^AMD{T^) + " " " • (5) 

Expectation values of a given tensor operator T^'^(rank A;) for the total-angular- momentum 
projected states P//^^^$^m£)(Z) and P^fj^^^ amd^) are calculated as follows, 

^pi},KMT^l\PtkK,^^) (6) 

= '^-^^^{J2MM\hM^) Y.{J2Kku\J,K^) J dQDi%{Q){<l>,\Tl;R{Q)\<^,). (7) 

The three-dimensional integral can be evaluated numerically by taking a finite number of 
mesh points of the Euler angle fl={a, 



B. Energy variation 

We make variational calculations for a trial wave function to find the state which mini- 
mizes the energy of the system; 

($1$) ■ ^ ^ 

In the AMD framework, the energy variation is performed by a method of frictional cooling, 
one of the imaginary time methods. Regarding the frictional cooling method, the reader is 
referred to papers |]2l|j2^ . The time development of the parameters Z of a wave function 
$(Z) is simulated by the frictional cooling equations, 

dXnk ^l d ($(Z)|g|$(Z)) . 1 . , 1 

"^^^^^^^^^ (^(Z)I^(Z)) ' = '^''""^ 

with arbitrary real numbers A and /i < 0. It is easily proved that the energy of the system 
decreases with each time step due to the frictional term://. After sufficient cooling iterations, 
the parameters for the minimum-energy state are obtained. 
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C. Wave function for states 

In order to obtain the wave function for a state, we make the energy variation after 
a spin-parity projection(VAP) for an AMD wave function by using the frictional coohng 
method explained above. That is to say we perform the energy variation for the trial function 
$ = P//^,$AAfD(Z), the spin-parity eigenstate projected from an AMD wave function. 
First we make variational calculations after only the parity projection but before the spin 
projection(VBP) to prepare an initial trial wave function ^AMDC^mu)- After obtaining an 
initial wave function in VBP calculations we evaluate the expectation values of Hamiltonian 
for the spin-parity-projected states by choosing the body-fixed 3-axis for the Q rotation to be 
the approximate principal z-axis on the intrinsic deformation. Then we find an appropriate 
K' quantum that gives the minimum diagonal energy of the spin-parity eigenstate 

where K' = (J3). For each spin parity J^, we start VAP calculations for the normalized 
energy expectation value (Pi^=^,$AA/D(Z)|i7|P//^,<l>AMD(Z))/(P//^,(Z)|P//^,(Z)) with the 
adopted K' quantum from the initial state. In general, the direction of the approximately 
principal 2;-axis is automatically determined in the energy variation because the shape of 
the intrinsic system can vary freely. The approximately principal axis can deviate from the 
3-axis in the VAP procedure with a given K' = (J3). That is to say that the optimum 
state P^i^AMD obtained after the variation may contain so-called i^'-mixing {K = {Jz)) 
components. However, the deviation of z-axis from the 3-axis is found to be small in many 
cases. It means that the obtained states do not contain i^'-mixing components so much, and 
K' well corresponds to the K quantum. 

Concerning the states in the lowest band with each parity, we can obtain appropriate 
initial wave functions by simple VBP calculations as mentioned above. For the highly 
excited states, in order to obtain initial wave functions and appropriate K' quanta for VAP 
calculations we make VBP calculations with a constraint on AMD wave function. The details 
of the AMD calculation with constraints have been described in Refs. In the present 



calculations, we adopt a constraint as the expectation value of the total-oscillator quanta to 
equal with a given number. After choosing a corresponding K' quantum, we perform VAP 
calculations from the initial wave function. When we obtain other local minimum states 
than the obtained states, we consider them as the states in higher rotational bands. 



D. Diagonalization 

After the VAP calculations for states, the optimum intrinsic states ^\md^ ^"amd^ ' ' 
and ^^MD obtained. Here m indicates the number of the calculated levels, we consider 
that the obtained wave functions approximately represent the intrinsic wave functions of the 

states. We determine the final wave functions by superposing the obtained AMD wave 
functions. That is to say that we determine the coefficients c, c', - ■ ■ in Eq. ^for each state 
by diagonalizing the Hamiltonian matrix {Plf^'^AMD\H\Pin<''^AMD) norm matrix 

{PMk'^AMD\PMk''^AMD) simultaneously with regard to {K\ K") and («, j). In comparison 
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with the experimental data such as the energy levels and the strength E2 transitions, the 
theoretical values are calculated with the final states after diagonalization. 



III. INTERACTIONS 



The adopted interaction for the central force is the case 3 of MVl force [2^, which 
contains a zero-range three-body teicm:V^^^ in addition to the two-body interaction:V^*^^\ 



Vdd = V^^^ + 

= -m + bP,- hP, - mP^P,) {vAexp [- 



+ Vr exp 



(12) 
(13) 
(14) 



where P^ and Pr stand for the spin and isospin exchange operators, respectively. As for the 



two-body spin-orbit force: Vl5, we use the G3RS force [27| as follows 



Vls = E^<J {uiexp i-Kjr^) + unexp {-Kjjr^)} (l+^^i^+fill . (s, + s,), (15) 

(16) 



The Coulomb interaction: Vc is approximated by a sum of seven Gaussians. The total inter- 
action V is a sum of these interactions: V = Vdd + Vls + Vc- 



IV. RESULTS 

The structures of the excited states of ^^Be are studied based on the VAP calculations 
in the framework of AMD. In this section we present the theoretical results concerning the 
energy levels, the E2 transitions, and the (3 transitions, which should be directly compared 
with the experimental data. More detailed discussions of the intrinsic structures are given 
in the next section. 

We adopt two sets of the interaction parameters. One parameter set (1) is m = 0.65, 
b = h = for Majorana, Bartlett and Heisenberg terms in the central force and uj = —uu = 
3700 MeV for the strength of the spin-orbit forces, which have been used in the previous 
study on the excited states of ^°Be. We also try another set (2) with weaker spin-orbit forces 
as uj = —uji = 2500 MeV. Other parameters in the case (2) are same as those in the case 
(1). The width parameter u is chosen to be 0.18 fm~^ which gives the minimum energy of 
^^Be in a VBP calculation. 

In VBP calculations we know that the lowest positive parity band K'" = 1/2+consists of 
1/2+, 3/2+,5/2+,7/2+ and 9/2+ states, and the lowest negative parity band K"" = 1/2" con- 
sists of 1/2^, 3/2^, and 5/2~ states. Therefore the jf" states in the lowest bands are obtained 
by VAP calculations for PI^,(^amd with the corresponding (J+, K') values as (1/2+, 1/2), 
(3/2+, 1/2), (5/2+, 1/2), (7/2+, 1/2), (9/2+, 1/2), (1/2", 1/2), (3/2", 1/2), (5/2~,l/2). We 
calculate the higher excited states in the second negative parity band by VAP calculations 
with (J+,ir')= (7/2-, 3/2), (9/2", 3/2), (11/2", 3/2), (13/2", 3/2), ■ ■ -. These states are 
considered to belong to a band with K'^ = 3/2~ . After obtaining the intrinsic states in 
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this second negative parity band for J > 7/2, the excited = 3/22 and 5/22 states in 
K'" = 3/2~ are found as local minimums with VAP calculations by starting from the intrin- 
sic states obtained for the higher spin states. In order to find the other excited 3/2~ state, 
we make VAP calculations for the spin-parity projected AMD wave function with the fixed 
intrinsic spin directions as 3 spin-down protons, one spin-up proton, 4 spin-down and 4 spin- 
up neutrons. The obtained wave function for the 3/2" state is dominated by a component 
of total-intrinsic spin Sp = 1 for protons. We superpose the wave functions to calculate the 
final wave functions by diagonalizing Hamiltonian and norm matrixes simultaneously. 

We should notice that all possible excited states are not exhausted in the present calcu- 
lations. In this work we perform VAP calculations basically for the rotational states which 
can be known from VBP calculations with or without constraint. For the higher excited 
states we must extend the VAP calculations with orthogonal condition to the lower states 



by superpositions as the previous studies [|T0|,^. Above the states shown in the present 



results, there should exist other excited states which may be obtained with the extended 
VAP calculations. 



A. Energies 

The theoretical binding energies of ^^Be are 58.2 MeV in the case (1) and 54.4 MeV in 
the case (2), both of which underestimate the experimental value 65.48 MeV. The binding 
energy can be reproduced by choosing a interaction parameter such as m = 0.60 of the 
Majorana term. However unfortunately it is difficult to reproduce all the features of nuclear 
structures such as binding energies, energy levels, radii, deformations and so on with one set 
of interaction parameters. Since we study the excited states taking care about the excitation 
energies and the intrinsic structures, we adopt the case(l) interaction which reproduces well 
the features of the excited states of ^°Be [jlO[ except for the binding energy. With this 



interaction the spin-parity 1/2+ of the ground state of ^^Be can be described in the present 
calculations. We also use another interaction case (2) with the weaker spin-orbit forces to 
be compared. The improvement of effective interactions is one of the important problems 
in nuclear studies. 

Here we comment on the stability of AMD wave functions above threshold energies. Since 
the single particle wave functions are written by Gaussians in AMD calculations, the relative 
motion between particles in a system is restricted by a Gaussian or a linear combination of 
Gaussians. Because of the limitation of the model space, continuum states nor out going 
waves can not be represented in the present model. Even if the energy of a nucleus is above 
the threshold energies of particle decays, the particles can not necessarily go away in the 
present framework. In this sense, the system is treated in the bound state approximation. 
The widths for the particle decays should be discussed carefully by other frameworks such 
as the method with reduced width amplitude or the complex scaling method beyond the 
present AMD framework. 

The energy levels of ^^Be are shown in Fig.^. There are many low-lying levels in the 
experimental data. The abnormal spin parity 1/2+ of the ground state has been known, 
although the normal spin parity of ^^Be is l/2~ in the simple shell-model picture. The 
calculations with the case(l) interaction reproduce the parity inversion between 1/2+ and 
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FIG. 1. The excitation energies of the excited states of ^-""Be. Theoretical results with the 
casc(l) and case(2) interactions are compared with the experimental data quoted from 'Table of 
Isotopes'. 
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l/2~ states. In the theoretical results, the rotational bands K'^ = 1/2+ and K'^ = l/2~ 
start from the band head 1/2+ and l/2~ states, respectively. Also in the results with the 
case(2) interaction, the 1/2+ state is lowest in the K'^ = 1/2+ band, however, it is slightly 
higher than the negative parity 1/2^ state by 1.4 MeV. As a result, the abnormal parity 
of the ground state can not be reproduced by the case (2) interaction which the spin-orbit 
forces are weaker than those in the case (1). We should not conclude the calculations with 
the case(l) interaction are better than those in the case (2), because the neutron halo effect 
on the energy gain of 1/2+ state is not taken into account in the present calculations. The 
details of the reproduction of parity inversion with the case(l) interaction are described in 
the later section. 

In both results in the case (1) and the case (2), there exist some rotational bands in 
the low energy region. By classifying the calculated excited states we can obtain three 
rotational bands K = 1/2+ , K = l/2~ and K = 3/2~ which are dominated by Ihu, Ohu 
and 2huj neutron configurations, respectively. In Fig.^, we show the excitation energies of 
the rotational bands as a function of total spin J(J + 1). We find a new eccentric band 
K = 3/2~ which starts from the second 3/2~ state at about 4 MeV. The excited states 
in the K = 3/2~ band are dominated by 2hijj excited configurations with 2 particles and 
3 holes in neutron shells. The K = 3/2~ band is constructed by a well-deformed intrinsic 
state with a developed clustering structure. The rotational band indicates large moment of 
inertia and it reaches the high spin states about 20 MeV. The highest spin is 15/2~ in the 
case (1) and 13/2^ in the case (2). In VAP calculations with further high spins, we can not 
obtain any stable states since an a particle escapes far away in the variational calculations. 
It means that the energy of the relative motion between clusters is beyond the classical 
barrier due to the coulomb and centrifugal forces. Since the barrier hight is very sensitive 
to the binding energy, we examine the highest spin of the K = 3/2~ band with another set 
of the interaction parameters taking care of the binding energy of ^^Be. If we change the 
Majorana exchange term of case (1) to m = 0.60, we can make the system bound as deeply 
as the experimental binding energy. Then a 17/2~ state is obtained as the highest spin 
state. W. Von Oertzen et al. suggested some candidates for the states belonging to this new 
negative parity band in the experimental data observed in the 2 neutron transfer reactions, 
^Be(^^C, ^^C)^^Be 0. Although they suggested a possibility of 19/2" state, the present 
results are negative to such a high spin state as 19/2" in the rotational K = 3/2~ band, 
because 17/2~ is the highest spin made from the 2p-3h configurations. It should be noticed 
that the present work is the first microscopic calculation which predicts the = 3/2~ band 
with the well-developed clustering structure. Since ^'^Be is a loosely bound system, further 
studies taking account of widths of the excited states are necessary to determine the band 
terminal. 

We find a excited 3/2~ state at about 10 MeV (Fig 0). In this state, the proton structure 
is quite strange comparing with the other excited states of ^^Be. In most of the states of 
^^Be, we find 2a cores as well as other Be isotopes: ^Be, ^Be and ^°Be. However, in the 3/2~ 
state at 10 MeV, only one a-cluster is formed. The other 2 protons do not form a a-cluster 
but couple to be totally S = 1 with aligned intrinsic spins. It is analogue with the structure 
of the first 1+ state with the unnatural spin parity in ^^C. We consider that this 3/2~ state 
must be the newly measured state at 8.04 MeV to which the /3-decay transition from 
^^Li is strong. Although its excitation energy is overestimated in the present calculations, 
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it is easily improved by changing parameters b and h of the Bartlett and Heisenberg terms. 
For example the parameters b = —0.2, h = 0.4, m = 0.41 gives 2 MeV lower excitation 
energy of the 3/2~ state than the one with b = 0.0, h = 0.0, m = 0.65. Here we fit the 
Majorana parameter m as to give the same a-a interaction as the case (1). The change of 
parameters gives no significant effect on the excitation energies of the other states with 2-a 
clusters. The results indicate that the Bartlett and Heisenberg terms should be taken into 
consideration in the detailed study of the energy levels though they have been often omitted 
in the traditional works on stable nuclei. 



B. Transition strength 

Data concerning the /?-decay strength are very useful to investigate structures of excited 
states. There are many experimental data concerning the f3~ and (3~^ decays into the excited 
states of ^^Be. The strength of the f3~ decays from ^^Li has been measured recently [Q. 
For the f3~^ decays, the strength of Gamov-Teller(GT) transitions has been deduced from 
the charge exchange reactions ^^B(t,^He)^^Be In the GT transitions from ^^Li(3/2~) 
and ^^B(3/2~), the allowed daughter states are l/2~, 3/2^ or 5/2^ states. In Table H the 
experimental log(/t) values are presented comparing with the theoretical results. In order 
to calculate the log(/t) values of GT transistions into the excited states of ^^Be we prepare 
the parent states ^^Li(3/2~) and ^^B(3/2~) by VAP calculations in the same framework. 

The experimental log(/t) values are reproduced well by the theoretical calculations. It is 
easily understood that the decays into the excited states in K = 3/2~ band are weak because 
these states have the ell-developed clustering structures with dominated 2huj components 
which make the overlap of a GT operator with the parent state of ^^B to be small. On the 
other hands, we consider that three levels at 0.32 MeV, 2.69 MeV and 3.96 MeV measured 
in the Z?"*" transitions correspond to the states 1/2 J", 3/2^" and 5/2^" in the K = l/2~ band 
because the experimental data indicate the significant strength as log(/it) < 5.0, which 
consistent with the theoretical results. 

What is important in the (3~ decays from ^^Li is that the strength is very sensitive to the 
breaking of 2 a-cluster cores in ^^Be. If the daughter state of ^^Be poses 2 ideal a-cluster cores 
with (Os)^ configurations, the GT transitions from ^^Li are forbidden completely because 
of Pauli principle. In another word, the strength of GT transitions from ^^Li indicates the 
degree of the a-cluster breaking in the daughter states of ^^Be. From this point of view, one 
of the reasons for the weak /5~ transitions as log(/t)=5.67 to the lowest 1/2^ state at 0.32 
MeV is the 2a-cluster structure in ^^Be(l/2~). Another reason for the weak f3~ transitions 
has been suggested to be the effect of the halo structure in ^^Li by T. Suzuki et al. pS 



The ground state of ^^Li is known to have the neutron halo structure which originates from 
s-orbits. The mixing of the s-orbits causes weak (3~ transitions to the normal states of ^^Be. 
Since the halo structure of ^^Li can not described in the present AMD wave function, the 
possible halo effects on the P~ decays are not included in the present results. If we suppose 
the mixing ratio of s-orbits in the ground ^^Li state as to be 50 %, the theoretical log(/t) 
values concerning the /3~ decays into l/2|f , 3/2^ and 5/2|f is expected to be increased by 
log(2)~ 0.3 due to the mixing. Even if we add 0.3 of the halo effect by hands to the present 
log(/t) values for 3/2{ and 5/2^^ states, the log(/t) values in case(l) are still smaller than 
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5.0, because of the sufficient cluster breaking. They are consistent with the experimentally 
measured rather small log/t values for the decays to ^^Be(2.7 MeV) and ^^Be(3.9 MeV). In 
another word, the experimental data indicates the significant breaking of 2 a-cluster cores 
in these states(2.7MeV and 3.9MeV) which we consider as the excited states 3/2^" and 5/2j". 

In the recent experiments of the measurements of l3~ decays, a new excited state at 8.04 
MeV with strong i3~ transitions has been discovered. The calculated log(/t) values for the 
3/2~ state of ^^Be at about 10 MeV well corresponds to this newly observed state at 8.04 
MeV. The transition is strong because the 2a-cluster structure disappears completely in this 
state. 

As shown in Table |, the experimental data concerning the strength of /3~ decays have 
been systematically reproduced in the present calculations. Considering the increases of 
log ft due to the halo effect of ^^Li, the theoretical values of log ft for j3~ calculated with 
the case (1) interaction well agree to the experimental data. The reason for the good 
reproduction is because the significant breaking of 2-a clusters in ^^Be is described in the 
present calculations. It is the same reason as the previous studies of ^^C in which the 
experimental data concerning the strength of l3~ and (3^ decays have been well reproduced. 
The quantitative discussion of the 2a-core breaking will be given in the next section. 

In Table we show the theoretical B{E2) values. It has been proved that AMD 
calculations reproduce well B{E2) values of light neutron-rich nuclei. As shown in the 
previous studies P, |22| , |23| , the experimental Q-moments and B{E2) values of various nuclei 
have been reproduced well by using the bare charges in the AMD framework because of 
the advantage of the AMD method which can directly describe proton deformations. Also 
in the present results, the features of proton-matter deformations in the intrinsic states 
are reflected in the theoretical B{E2) values. In the lowest negative parity band = 
l/2~, the intrinsic system deforms because of 2 a-cluster cores although the development of 
clustering is smallest comparing with the other bands. B{E2) values are larger in the ground 
K'" = 1/2+ band which has the developed clustering structure with a large deformation. The 
B{E2) values are enhanced between the states in the second negative parity band K'" = 3/2~ 
with an extremely large deformation due to the clustering development. 

The El transition strength from 1/2^" to 1/2^ in ^^Be is known to be very large comparing 
with those in other light nuclei. The theoretical value of B{E1; l/2i l/2]'")=0.02 e-fm 
is much smaller than the experimentally observed strength 0.116 ± 0.011 e-fm. The halo 
structure may give an important effect on the El strength. 

V. DISCUSSIONS 

We investigate the structures of ^^Be focusing on the clustering aspects. It is found 
that two a-cluster cores are formed in most of the states of ^^Be as well as in other Be 
isotopes:^Be, ^Be, ^°Be. In the present calculations, we found three rotational bands with 
2-a cores in ^^Be. In this section, we argue the development of clustering and investigate 
the roles of valence neutrons in the clustering states. We analyze the breaking of cluster 
cores and also describe the feature of a non-clustering state. The problem of parity inversion 
between positive and negative parity states is also discussed. 
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TABLE I. Log(ft) values of the P transitions. The theoretical values are obtained from the 
Gamow- Teller transition strength. The experimental data of /? decays from ^^Li are taken from 
Ref. Q. Log ft values concerning the Gamow- Teller transition from ^^B are deduced from the 
charge exchange reactions 
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TABLE II. E2 strength of ^^Be. The theoretical results of AMD with the interactions case(l) 
are listed. 
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A. Clustering structure 

We investigate the intrinsic structures of the excited states focusing on the clustering 
aspects. Although the states PIik^^md superposed for all the obtained wave functions 
^\md so as to diagonalize a Hamiltonian matrix, the J+ state after the diagonalization is 
found to be dominated by the AMD wave function PitK^^MD which is obtained in a VAP 
calculation for the given spin and parity J+. Therefore we consider an obtained AMD wave 
function ^amd with a VAP calculation as the intrinsic wave function of the corresponding 
J+ state. 

In the excited states, three rotational bands: K = 1/2+ , K = 1/2- and K = 3/2" 
are recognized, because the intrinsic structures of the states in each band are similar. The 
density distributions of the intrinsic states of the band head states: 1/2^^, 1/2^ 3/22 are 
shown in Fig.^. The neutron structure of an intrinsic state is very different from the ones of 
the other rotational bands. The density distributions of neutrons are presented in the right 
column in Fig.^. As for the neutron deformation, the 1/2" state has an oblate deformation. 
The oblate shape is natural for the normal parity state in a system with the neutron number 
N = 7. In the positive parity 1/2+ state, the neutron density deforms prolately. The prolate 
deformation of neutron density is extremely enhanced in the 3/22 state. The detailed 
discussions of neutron structures is given later by analyzing the single-particle behavior of 
neutrons. 

In the proton density shown in the middle column in Fig.|^, we can see dumbbell shapes 
due to 2 pairs of protons. Roughly speaking, it indicates that 2-a cores are formed in all 
the intrinsic states of these three rotational bands. The spatial development of clusters is 
small in the lowest negative parity 1/2]^ state which has a main component of a normal 
Ohuj configuration. In the ground 1/2+ state with the abnormal parity, the 2a clustering 
structure develops. It is interesting that the most remarkable clustering structure is seen 
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FIG. 3. The density distributions of the intrinsic structures of 1/2^", 1/2+ , 3/22 states 
calculated with the case(l) interaction. The intrinsic system is projected on to a plain which 
contains the longitudinal axis of the intrinsic states. The density is integrated along a transverse 
axis perpendicular to the plain. The density distributions for matter, proton and neutron are 
presented in the left, middle and right, respectively. 

in the 3/22 state which belongs to the K — 3/2~ band. In this state, the relative distance 
between 2-a cores is largest among these 3 rotational bands. It seems that the 2a-clustering 
structure develops following the prolate deformation of neutron density. According to the 
development of 2a clustering, the proton deformation is larger in the = 1/2+ band than 
in the K — 1/2" band, and largest in the K — 3/2~ band. The difference of the proton 
structures are reflected in the theoretical B{E2) values as mentioned in the previous section. 
Namely, the theoretical values concerning the strength B{E2) are large in the transitions 
between the states in the rotational bands with the large proton deformations. As mentioned 
above, the K = 3/2" band is the largely deformed state with the well development clustering 
structure. However, the deformation decreases and the clustering weakens at the band 
terminal state with the increase of spin near the highest spin J = 15/2~ in case(l). 

B. Behavior of valence neutrons 

As mentioned in the previous subsection, the deformation of neutron density changes 
drastically between the rotational bands. We inspect the single-particle wave functions of 
the intrinsic states to study the behavior of valence neutrons. The single-particle wave 
functions of an intrinsic state are obtained by transforming from the original Gaussian 
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functions to the orthonormal single-particle wave functions under the diagonal condition of 
the single-particle Hamiltonian matrix. The definition of single-particle wave functions and 
energies of an AMD wave function is described in Refs. PJl^ • 

In Be isotopes, it has been known that the molecular orbits around la play important 
roles. The molecular orbits in Be isotopes have been suggested in ^Be with 2a + n cluster 
models [0. W. Von Oertzen et al. proposed a picture of Be dimers to understand 



the excited states of neutron-rich Be isotopes. They predicted ex- and 7r-orbits which are 
made from linear combinations of p-orbits surrounding 2-a clusters (see Fig.^. Itagaki 
et al. have described the structures of ^°Be and ^^Be with a extended cluster model by 
assuming a cores and the molecular orbits. In the study of Be isotopes with AMD methods 
i8|j9|j2^, the structures with 2 as and valence neutrons in neutron-rich Be isotopes have 
been microscopically confirmed without assuming the existence of any clusters nor molecular 
orbits. In the present work we investigate the behavior of the valence neutrons focusing on 
the molecular orbits. 

By analyzing the ratio of positive and negative parity components in each single-particle 
wave function, it is found that all the neutron orbits are approximately parity-eigen states, 
which roughly correspond to Is-, Ip-, and 2s\d orbits. The lowest normal-parity l/2~ state 
comes from a O^u; configuration because 2 neutrons occupy s-like orbits and the other 5 
neutrons are in p-like orbits. On the other hand, in the ground 1/2"'" state with the ab- 
normal parity, the last neutron occupies a srf-like orbit, which means that the 1/2+ state 
is dominated by a XTiio configuration. An interesting feature is found in the highest neu- 
tron wave function which well corresponds to the S(i-like orbit. The density distribution 
of neutrons are presented in Fig.^(a). It has nodes along the longitudinal direction of 2-a 
clusters and contains nearly 90% of a positive-parity eigen state. It well corresponds to the 
so-called molecular cr-orbit explained in the schematic figure ^(a). On the other hand, the 
features of 2 neutron orbits energetically below the cx-orbit are shown in Q^b). The domi- 
nating negative-parity component of these orbits indicates that the orbits correspond to the 
p-like orbits. As shown in Fig.^(b), the orbits seem to be similar to the molecular vr-orbits. 

One of the reasons for parity inversion of ^^Be is considered to be because of the molecu- 
lar cr-orbit which originates from sd-orbits. Since Be isotopes prefer to prolate deformations 
because of 2-a clusters, the cx-orbit can easily gain its kinetic energy in the clustering de- 
veloped system. In another word, the 2a-clustering structure is one of advantages for the 
cx-orbit in Be isotopes. From this point of view, it is easy to understand the reason for 
the clustering development in the 1/2^ state in terms of the energy gain of the cx-orbit in 
a clustering system. Figures |^(a),(b),(c) show the single-particle energies in the intrinsic 
states of l/2j", 1/21" and 3/22 with the case (1) interaction, respectively. The highest orbit 
in Fig.p(b) corresponds to the cx-orbit in 1/25*" state. The energy of the intruder cx-orbit is 
lower than the highest vr-orbit in 1/25". ^ the states: 1/2^, 1/2^ and 3/22 , the 2-a cores 
consist of 4 neutrons and 4 protons in the Is-like orbits and the lower p-like orbits seen in 
Fig.^. The behavior of the single-particle wave functions of valence neutrons in the present 
results is almost consistent with the previous AMD studies 

One of new discoveries in the present work is that the possible negative-parity band 
K = 3/2~ has a largely deformed shape with a remarkable clustering structure. With 
the same analysis of the single-particle wave functions as mentioned above, we find that 2 
neutrons in the 3/22 state occupy so-called cx-orbits. In another word, the well-developed 
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FIG. 4. Sketches for the molecular orbits, (a) u-orbits and (b)7r-orbits surrounding 2as. The 
molecular orbits are understood by linear combinations of p-shell orbits around the a clusters. 



o 

X (10 fm) X (10 fm) 

FIG. 5. The density distribution of single-neutron wave functions in the intrinsic system of 
1/2]'' state calculated with case (1). Figure (a) shows the density for the highest neutron orbit 
with 90% positive parity component, while figure (b) is for the second neutron orbit which contains 
90% of a negative parity component. 




clustering is caused by these 2 neutrons in the a-orbits so as to gain their kinetic energies. 
It is surprising that such a 2 fiio state appears in the low energy region. The clustering 
structure in the K = 3/2^ band is the very molecular structure which has been predicted 
by W. Von Oertzen 0. This is the first microscopic calculation which shows that such an 
exotic structure should exist in the low energy region. 

In the single-particle energies shown in Fig.p(c), we find another interesting thing that 
the energies of the cx-orbits in 3/2^ state become lower than that of an odd neutron 7r-orbit. 
One of reasons for the higher energy of the odd neutron vr-orbit than the cr-orbits is a lack of 
the pairing interaction because of the absence of the other neutron in the vr-orbit. Another 
reason is the degeneration between cr-orbits and vr-orbits at the large inter-cluster distance. 
With the increase of relative distance between 2-a clusters, the cr-orbits gain their energy, 
while the vr-orbits loose their potential energies. As the relative distance increases, the cr- 
orbits cross the higher vr-orbits and come down near the lower vr-orbits. It is helpful to see 



the neutron level scheme of two center shell model |^1| to understand the inversion of vr- 
and the cr-orbits with the increase of inter-cluster distance. 

In the single-particle energies shown in Fig.|^, the small level spacing between the lower 
Is-orbits and the lower Ip-orbits indicates the development of 2 a clusters in 1/2+ and 3/2^ 
states. In fact, the level spacing is smallest in 3/22 state. 

As mentioned above, we see the drastic change of neutron structures between three 
rotational bands. It is found that the origins of drastic change are the neutrons in cr-orbits 
and vr-orbits surrounding 2-a cores. The intrinsic states are characterized by the number 
of valence neutrons in the cr-orbit. That is to say that, no neutron, one neutron and two 
neutrons occupy the cr-orbits in the K = 1/2' K = 1/2+ K = 3/2^ bands, respectively. 
The deformation with clustering structure is enhanced by the increase of the neutrons which 
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FIG. 6. Single-particle energies in the intrinsic system of (a)l/2^ , (b)l/2]'", (c)3/22 with the 
case (1) interaction. Energies for protons (neutrons) are presented in left(right) in each figure. 
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occupy the cx-orbits due to the energy gain of a-orbits in the developed clustering structure. 



C. Breaking of 2a cluster cores 

Although the structures with 2 a-cluster cores are seen in most of the states of ^^Be, 
the a clusters are not the ideal a clusters with simple (Os)^ configurations. The a clusters 
are broken slightly because of the spin-orbit force. The components of cluster breaking 
allows the /? decays from ^^Li into ^^Be. This is a similar situation as seen in the Sa-cluster 
structures of ^^C to which Gamow- Teller transitions from ^^B and ^^N are not weak. 

In the calculated results, we can know the degrees of breaking from the ideal 2-a clusters 
by estimating non-zero total-intrinsic spin(S'p 7^ 0) components for protons. If we neglect 
Sp>2 components and assume that only S'p = and Sp = 1 components are contained in a 
state, the value 0.5 x {S^) indicates directly the mixing ratio of the non-(Os)'^ configurations, 
which signifies the degrees of cluster breaking in the state. In table |T|, the expectation values 
of the squared total-intrinsic spin of protons are listed. Generally speaking, the clustering 
breaking is larger in the results in the case (1) than in the case (2) because the stronger 
spin-orbit forces in the case (1) give larger effects on the dissociation of 2 a-cluster cores. 
The states l/2j~ , 3/2j~ and 5/2j" in = l/2~ contain significant components of the cluster 
breaking in the ratios of 7%, 10% and 20%, respectively. The mixing of cluster breaking 
allows the (3 decays into these states from ^^Li, which consistent with the experimental 
measurements (Table.|). The smaller mixing in 1/2^ than the other 3/2^ , 5/2^ states 
is reflected on the relatively weak (3 decays to l/2|f as mentioned in the previous section 
on the /3-decay strength. On the other hand, the breaking of a clusters are very small in 
K = 1/2"*" and K = 3/2" which have the well-developed clustering structures. Therefore 
Gamow- Teller transitions to the excited states in K = 3/2~ band from ^^Li are predicted to 
be weak except for the /3 decays to the 3/22 state because this state slightly contains the 
cluster breaking due to the state mixing with the lower 3/2|f state. 

We discover an excited 3/23 state where the 2a-cluster structure is completely broken. 
In the density distribution of protons we can not recognize the dumbbell shape in this state 
(see Fig.|^). As seen in Table |T|, the main component is the proton intrinsic spin 5^ = 1 
state. As a result, the /3-decay transitions from ^^Li to the 3/23 state are strong comparing 
with the /3-decay strength to the other excited states of ^^Be. It is consistent with the strong 
f3 decays to the state at 8.04 MeV, which have been recently measured. 

Although the 2a-core structures are dominant in most of Be isotopes, it is interesting 
that the excited state with the completely breaking of one of 2 as appears in neutron-rich 
Be. Even in the excited states with 2-a cores, the components of the cluster breaking play 
an important role in the strength of (3 decays from ^^Li. 



D. parity inversion 

We consider details of parity inversion of the ground state in the case(l) calculations. 
Concerning the energy gain of the positive parity 1/2"*" state, the effects of the halo struc- 
ture, the clustering structure with a prolate deformation, the three-body force, the angular- 



momentum projection, and the paring have been discussed in the pioneering works |TB 
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TABLE III. The expectation values of the squared total intrinsic spin of protons in the 
excited states of ^^Be. 
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FIG. 7. The density distribution of the intrinsic wave function of 3/23 calculated with the 
case(l) interaction. 
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In the present calculations the effects of clustering structure are included as well as the 
pioneering AMD studies [0,^ • One of the reasons for parity inversion is the energy gain of 
the molecular cr-orbit in the system with developed la clustering. The details of clustering 
and molecular orbits have been already mentioned above. 

The three-body force was calculated approximately with a perturbative treatment in the 
pioneering study [0. In the present work, the three-body force is treated exactly as same 
way as the previous AMD study in Ref. 



In the previous AMD studies of Be |16,22], the angular-momentum projection was 



perturbatively treated. That is to say that the energy variation was performed after a parity 
projection but before an angular- momentum projection(VBP). In the present results, we 
perform energy variation exactly for the spin-parity projected wave functions(VAP). Then we 
find that the effect of angular-momentum projection on the parity inversion is quite different 
from the previous VHP results with perturbative treatments. In order to quantitatively 
discuss the effects, we estimate the energy difference between l/2~ and 1/2+ states: bE = 
Exj2+ ~Eii2-- The parity inversion is described by a negative value of 5E. The experimental 
value of 5E has been known to be —0.32 MeV. For simplicity we do parity-projected AMD 
calculations with fixed intrinsic spins by using the case(l) interaction. In table the values 
of 5E and energies (a) without angular-momentum projections, (b) the angular-momentum 
projections after variation within the perturbative treatments of VBP, and (c) results of 
the exact VAP calculations are listed. Before the angular-momentum projections, 5E is 0.6 
MeV. After the projection within the perturbative treatments of VBP, 5E becomes small 
as —1.1 MeV. It is because that the energy gain of the positive-parity state by the angular- 
momentum projection is 1.7 MeV larger than the gain of the negative-parity state. The origin 
is the well-developed clustering structure in the former state comparing with the latter state 
because a deformed system gains much energy due to a angular-momentum projection in 



general. The VBP results are consistent with the previous AMD studies in Refs. |]22| , p!6 
However when we exactly treat the angular-momentum projection with VAP calculations 
we find a new aspect in the effects on the energy gain. In the previous AMD studies, it 
is known that in many cases VAP results of low-lying states is not so much different from 
VBP results. It is true in the 1/2+ state. However it is not true in the 1/2^ state. In 
the VAP calculations, the negative parity state gains much energy as 8.6 MeV due to the 
angular-momentum projection because the deformation of l/2~ state grows more largely in 
the VAP calculations than in the VBP calculations so as to gain the energy of the angular- 
momentum projected state. As a result, 5E of (c) is 1.0 MeV as much as the result before 
projections (a). In fact, the clustering of ^^Be(l/2]^) is larger in the VAP results than as is 
expected in the VBP calculations. Therefore the effect of the angular-momentum projection 
is not positive for the parity inversion in ^^Be. This is one of new things in the present 
work which points that the perturbative treatments of angular-momentum projections do 
not always work well for precise descriptions. 

Comparing with the previous AMD studies of low-lying states of ^^Be ||T6|j2^ , another 
improvement in the present framework is treating the intrinsic spin functions as variational 
parameters. In the Tabled, the effect of the flexible intrinsic spins are seen in the difference 
between (c) the VAP calculations with the fixed intrinsic spins and (d) the VAP calcula- 
tions with the variational intrinsic spins. The treatment of flexible intrinsic spins gives an 
important effect on the parity inversion of ^^Be. It reduces 5E to a 1.1 MeV smaller value 



21 



TABLE IV. Energies of the positive and negative parity states of ^^Be. The parity projected 
AMD calculations with fixed intrinsic spins (a) before angular-momentum projections, (b) after 
angular-momentum projections of VBP within the perturbative treatment, and (c) results of the 
exact VAP calculations are listed. We also show (d) the present VAP calculations with variational 
intrinsic spins but no diagonalization of the basis. The adopted interactions are the case(l) force. 
The energy difference 5E between the positive and negative parity states is defined as £'1/2+ —E1/2- 
for the calculations after angular-momentum projections and E{+) — E{—) before the projection. 
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0.6 


-1.1 


1.0 
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in (d) than that in the calculations (c). The importance of the flexible intrinsic spins in 
P3/2 shell-closure states has been already argued in the studies of ^^C p3| by one of the au- 
thors. In the case of ^^Be, it is natural to consider that the sub-shell effect is one of the 
reasons for the parity inversion because the pi/2-orbit sould be raised relatively higher than 
2si/2-orbits. Unfortunately, the simple AMD method with the fixed intrinsic spins describe 
the sub-shell effects insufficiently. This is one of the advantages of the present framework 
which can describe well both aspects of the clustering and shell effects. 



Regarding the other reasons of parity inversion, Sagawa et al. and Dote et al. ||T6 
suggested the paring effects in the neutron p-shell. In the present framework, a part of the 
paring effects in p-shell should be included automatically by the spin-parity projections and 
the superpositions. 

As mentioned above, the clustering effects, the three-body forces, the P3/2 sub-shell effects 
are included in the present calculations. The paring effects are expected to be partially 
contained. The other effect which should be important for parity inversion is the neutron- 



halo effect According to the study in Ref. [|T^, the effect due to the neutron halo 



structure is estimated as 0.6 MeV reduction of 5E. However the halo effect is not taken into 
account in the present calculations. We think that this halo effect and the residual paring 
effect are contained effectively in the strength parameter uls = 3700 MeV of the spin-orbit 



force which is slightly stronger that the parameter u^s = 3000 MeV adopted in Ref. p3 



The halo structures should be important only for the loosely bound states with low-spin 
orbits such as s-orbits and j9-orbits at most. The present strong effective spin-orbit force 
in case(l) may artificially reduce too much the theoretical excitation energies of the other 
states in = 1/2+ band except for the 1/2"'" state. Before concluding the relative position 
of the energies between K = 1/2^ and K = 1/2+ bands we should more detailed researches 
by taking the remained effects such as the halo structure into consideration. 
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VI. SUMMARY 



We studied the structures of the ground and excited states of ^^Be with VAP calculations 
in the framework of antisymmetrized molecular dynamics. Various kinds of excited states 
with the clustering structures and also the non-clustering structures were discovered in the 
theoretical results. Wc predicted many excited states. Most of the states belong to three 
rotational bands: K = 1/2+ . K = 1/2^ and K = 3/2^, which are dominated by Ihu, 
Ohcu and 2hu! configurations, respectively. It should be pointed out that the formation of 2 
cc-cluster cores is seen in many excited states of ^^Be in the the present resutls in spite of no 
assumption of the existence of clusters. The interesting point is that an eccentric rotational 
band K = 3/2~ with the mostly developed clustering structure starts from the 3/2^ state 
at about 4 McV and reaches high spin states. 

The experimental data concerning the /S^-decay and /9~-decay strength were reproduced 
well. We also argued that the cluster breaking plays an important role to allow the (3 decays 
from ^^Li. The significant breaking of 2-a cores in the states in K — l/2~ band has been 
seen in quantitatively estimating the breaking of clusters. We discovered a non-clustering 
state at about 10 MeV. One of the characteristics of this state is the strong f3 decays from 
^^Li, which well corresponds to the new excited states at 8.04 MeV found in the /3~-decay 
measurements. 

By analyzing the single-particle wave functions in the intrinsic states, it was found that 
the molecular cr-orbits surrounding 2-a cores play important roles in the clustering structures 
of ^^Be. In the ground band K = 1/2"*", one neutron occupies the cr-orbits. In the newly 
predicted K = 3/2~ band is dominated by 2hu! configurations with two neutrons in the 
cr-orbit. When the surrounding neutrons occupy the cr-orbits, the clustering development 
is enhanced so as to gain their kinetic energy. In another word, one of the reasons for the 
parity inversion and the low-lying 2hu! states is the energy gain of the cr-orbits with the 
developed clustering structures. 

Concerning the mechanism of the parity inversion, we described the importance of molec- 
ular neutron orbits in the developed clustering structure and also mentioned about the P3/2 
sub-shell effect. Although the spin parity of the ground state is described with a set of 
interaction parameters case (1), we need more detailed researches taking the halo structures 
of the ground state into account. 
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